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We study the dynamics of the quantum phase transition of a ferromagnetic spin-1 Bose-Einstein 
condensate from the polar phase to the broken- axisymmetry phase by changing magnetic field, and 
find the spontaneous formation of spinor domain walls followed by the creation of polar-core spin 
vortices. We also find that the spin textures depend very sensitively on the initial noise distribution, 
and that an anisotropic and colored initial noise is needed to reproduce the Berkeley experiment 
[Sadler et a/., Nature 443, 312 (2006)]. The dynamics of vortex nucleation and the number of 
created vortices depend also on the manner in which the magnetic field is changed. We point out 
an analogy between the formation of spin vortices from domain walls in a spinor BEC and that of 
vortex- antivortex pairs from dark solitons in a scalar BEC. 

PACS numbers: 03.75.Mn, 03.75.Lm, 03.75.Kk, 67.57.Fg 



INTRODUCTION 



Topological defects have played a key role in under- 
standing the physics of scalar and spinor superfluids. 
The first topological defect observed in a gaseous Bose- 
Einstein condensate (BEC) was a quantized vortex cre- 
ated using the phase imprinting method [l|. Vortices 
have also been created using rotating potential 0] and 
by means of adiabatic spin rotation accompanied by a 
topological Berry phase 3]. In these methods, topo- 
logical defect formation is enforced by an external laser 
or magnetic field. On the other hand, topological de- 
fects can also be created spontaneously. The Mermin-Ho 
texture 3 of superfluid 3 He in a cylindrical container 
involves a coreless vortex. The Kibble- Zurek mecha- 
nism 0, |(| in a quenched superfluid or in the early uni- 
verse affords another intriguing example of spontaneous 
topological defect formation. Thermally nucleated vor- 
tices have recently been observed in a gaseous BEC, pre- 
senting yet another topological phase transition known 
as the Berezinskii-Kosterlitz-Thouless transition Q . 

Recently, the Berkeley group 8] observed the spon- 
taneous topological defect formation of spin in a spin-1 
87 Rb BEC. The atoms were prepared in the m = state 
(polar phase in Fig. in a strong magnetic field, where 
m is the magnetic quantum number. The magnetic field 
is then lowered below a certain critical value (broken- 
axisymmetry phase in Fig. [IJ, where the m = state 
becomes dynamically unstable and magnetization grows 
in a direction perpendicular to the magnetic field, and, 
consequently, the axisymmetry in the spin space is spon- 
taneously broken In the experiment, the magnetic 
field was rapidly decreased across the critical value (as 
shown by the arrow in Fig. — a process we refer to 
as "quench" — and the ensuing magnetization dynam- 
ics was observed by the spin-sensitive in situ measure- 
ment 10]. After the quench, it was observed that mag- 
netization grew to form complicated ferromagnetic do- 
mains, and that, remarkably, some snapshots of the spa- 
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FIG. 1: Phase diagram of an homogeneous spin-1 BEC with 
ci < 0. In the ferromagnetic phase, the ground state is the 
m — 1 state for p > and the m — — 1 state for p < 0, while 
in the polar phase, the ground state is the m — state. The 
ground-state wave function in the broken- axisymmetry phase 
possesses three nonzero components given in Eq. (0, and the 
spin vector tilts against the applied magnetic field. The arrow 
indicates the direction of the quench of the magnetic field. 



tial spin distribution revealed topological spin textures, 
known as polar-core spin vortices [llj . 

It has been predicted that a variety of spin textures, 
such as the staggered domain and helical textures, are 
created by the quench of the magnetic field The 
spontaneous nucleation of the polar-core vortex has also 
been predicted in Ref. 13]. The underlying physics in the 
spontaneous spin-texture formation is spin conservation, 
which prohibits uniform magnetization. 
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Motivated by the Berkeley experiment in the 
present paper we study the dynamics of a spin-1 BEC 
caused by the quench of the magnetic field from the 
polar to the broken- axisymmetry phase, and investigate 
the formation dynamics of the topological spin texture. 
We show that the spin vortices observed in Ref. are 
formed in two steps. First, spin domains separated by 
domain walls develop, breaking the axisymmetry in the 
spin space. Secondly, the domains transform into spin 
vortex- antivortex pairs. We find that the details of the 
dynamics depend on the initial seed in the m = ±1 states, 
which represents residual atoms, quantum fluctuations, 
and thermal noises, and identify an initial seed that re- 
produces the experimental results. We also study the 
quench-time dependence of the number of created spin 
vortices. 

This paper is organized as follows. Section ITU reviews 
the mean-field theory of the spin-1 BEC and some of 
its key properties relevant to later discussions. We also 
point out an analogy between defects in a spinor BEC 
and those in a scalar BEC in Sec. Ill Cl Section HTT1 stud- 
ies the magnetization dynamics of the trapped system. 
Sections 1111 CI 1111 Dl and 1111 El examine three different 
kinds of initial noise conditions and identify an appropri- 
ate one that captures the main features of the Berkeley 
experiment. Section All Fl investigates the dependence of 
the dynamics on the speed of the quench and on the fi- 
nal value of the magnetic field. Section ITlI Gl studies the 
behavior of the growth of the total magnetization, and 
Sec. IIVI concludes the paper. 



II. MEAN-FIELD THEORY FOR SPIN-1 BEC 

A. Formulation of the problem 

The zero-temperature mean-field energy of a spinor 
BEC confined in an optical trapping potential V is given 

by 



(1) 

where M is the mass of the atom and ip m is the con- 
densate wave function for atoms in magnetic sublevel m, 
satisfying J2 m I dr\ifj m \ 2 = N with N being the total 
number of atoms. We define the atomic number density 



and the spin density 
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where / = {fx-, fy, fz) is the vector of the spin-1 matrices. 
The interaction energy E- int for the spin-1 atom has the 



form QUI 
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with as being the s- wave scattering length for two col- 
liding atoms with total spin S. In the present paper, 
we consider the case of spin-1 87 Rb atoms, and take 
ao = 101. 8clb and <22 = 100. 4ae 16], where ae is the 
Bohr radius. The linear and quadratic Zeeman energy 
Eb under magnetic field B is given by 
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where hb is the Bohr magneton, is the hyperfine 
splitting energy, and 1/2 in the first term is the Lande 
g-factor. The mean-field dynamics of the system is thus 
described by the multicomponent Gross-Pitaevskii (GP) 
equation, 
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where q = (/^bB) 2 /(4E h f), F± = F x ±iF y , and the mag- 
netic field is assumed to be in the z direction. 

The linear Zeeman terms ^/hbB/2 in Eq. {Taj only ro- 
tate the spin around the z axis at the Larmor frequency. 
Going onto the rotating frame of reference by setting 

^±1 e TWBBt/(2h)^ ±i and F± _^ e ±in B Bt/(2h) F± ^ we 

find that the linear Zeeman terms can be eliminated. 



B. Ground state and excitation spectrum 

We briefly review the ground state and excitation spec- 
trum for the homogeneous case 0,EJ. For spin-1 87 Rb 
atoms, ci is negative, and the ground-state phase dia- 
gram is given as Fig. ^ where p = /hbB/2 + po with po 
being the Lagrange multiplier [l7|. Here, po is introduced 
to set F z at a prescribed value, since the total magnetiza- 
tion in the z direction is conserved. When the quadratic 
Zeeman energy q is larger than 2|ci|n, the polar phase, 
C01? ^Oj ip-i) = V^(0> 1 5 0), is the ground state. Between 
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the polar and ferromagnetic phases, there is a broke n- 
axisymmetry phase shown as the shaded region in Fig.^] 
The order parameter in this phase is given by 9] 
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where x±i are arbitrary phases of the m = ±1 states. 
We note that spin vector F in this phase tilts against the 
direction of the magnetic field, breaking the axisymmetry 
spontaneously 9]. The phase boundaries are given by 
q = \p\ and 2\c\\nq — q 2 +p 2 = 0, across which the system 
undergoes the second-order phase transition. 

In Sec. IIIII in an attempt to study the Berkeley ex- 
periment we will examine the magnetization dynamics 
from the initial m = state. In order to understand the 
dynamics qualitatively, we consider the Bogoliubov exci- 
tations from the m = state for the homogeneous case. 
Solving the Bogoliubov-de Genne equation, we obtain the 
excitation spectrum as 0, 



Eo = y/e k (e k + 2c n), 



(9) 



E ±1 = T^bB + y/(e k + q)(e k + q + 2cm), (10) 

where e k — h 2 k 2 /(2M) with k being the wave number. 
The mode given in Eq. (0) involves only the m = state, 
and can be regarded as a phonon mode, whose excitation 
energy is always real and positive. The two modes given 
in Eq. ((10)) are magnon modes, which transfer the atoms 
from the m = state to the m = ±1 states. The excita- 
tion energies E±\ are complex for — q < e k < —2c\n — q. 
Therefore, when c\ is negative and q < 2|ci|n, the m = 
state is dynamically unstable against long- wavelength ex- 
citations of the magnon modes. For q < |ci|n, the most 
unstable wavelength is 
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and the corresponding imaginary part of E±\ is |ci|n, 
giving a characteristic time scale for the dynamical in- 
stability as 

h (12) 



7~mu — 



\ci\n' 



For a larger magnetic field satisfying \c\\n < q < 2|ci|n, 
the most unstable wave number is k = 0, and there- 
fore the system tends to magnetize uniformly. The cor- 
responding time scale is given by Ti/[q\q + 2cin|] 1 / 2 . 

The above Bogoliubov analysis is valid only when the 
system is homogeneous and deviations from the initial 
state are small. To investigate regions beyond such re- 
strictions, we will solve full GP equation (0 numerically 
in Sec. IIIII and analyze the detailed magnetization dy- 
namics of the trapped system. 



C. Topological defects in spinor BECs 

Quantized vortices and dark solitons in scalar BECs 
are topological defects, in which the density vanishes due 
to the topological constraint on the phase of the wave 
function. We relate these topological defects in scalar 
BECs to spin vortices and domain walls in spinor BECs, 
in which local magnetization vanishes. 

We first note that throughout the spin-exchange dy- 
namics the total density n = | ^_ 1 1 2 + 1 | 2 + 1 V^i 1 2 remains 
almost constant because Co ^> \c\\. Since the transverse 
magnetization develops from the m = state, we also as- 
sume that |^-i | = |^i I an d phases Xm of ip m are related 
to each other by x±i = Xo ± where a is an ar- 

bitrary function. Substituting these relations in the GP 
equation 0, we can eliminate ifj-i and ^o, obtaining 
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where we drop the trapping potential and the linear Zee- 
man term for simplicity. Equation (|13|) is the single- 
component GP equation with a repulsive interaction 
(-4ci > 0). 



Polar- core vortices 



The general form of a spin- vortex state is given by 



V-m(r) =e i ^e ic ^f m (r ± ), 



(14) 



where <j) is the azimuthal angle, 7 m 's are global phases, 
c m 's are integers, and / m 's are real functions of r±_ = 
(x 2 + y 2 ) 1 / 2 . In Eq. {T2J), we assume that the system 
is uniform in the z direction and that the vortex core 
is located on the z axis. For example, the Mermin-Ho 
texture has topological charges c\ = 0, Co = 1, and c_i = 
2. 

The solution to Eq. ((13)) with a singly-quantized vortex 
is given by 



/o(r±), 



fi(r±), 



(15) 



where 7 is an arbitrary phase and f§(r±) + 2f 2 (r±) = n. 
Since the m = ±1 states have topological defects, fi(r±) 
vanishes at ri = 0, and the core is then occupied by 
the m = (polar) state. The state given by Eq. ([T5| 
is therefore referred to as a polar-core vortex. Equation 
(fill)) indicates that f\ is proportional to ri/( s near r± = 
0, where 
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and hence the size of the vortex core is characterized by 

6- 
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The spin density of Eq. llTnj) is given by 

F + = V2(r^o + tl>W-i) = 2V2e^- i y 1 (r ± )f (r ± ), 

(17) 

F z = 0. Equation JT7J) shows that the spin vector cir- 
culates around the z axis with phase winding =pl and 
vanishes at r± = 0. In fact, the spin current 
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for Eq. ([T5|) is calculated to be Jf- = J s y • = and 
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where is a unit vector in the azimuthal direction. On 
the other hand, state JTSJ) has no mass current, 
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2Mi 



E^Wm-(V^ m ], (20) 



and has no orbital angular momentum. 



2. Domain walls 



Equation ([T3| has a dark-soliton solution as 



l2\ Cl \n- 
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where \i = cqu— \c\ \n+q/2 and we assume that the planar 
dark soliton is located at z = 0. This solution describes a 
domain wall at z = 0, since the spin vector F vanishes at 
z = and asymptotically approaches a constant vector 
with opposite directions for z — > ±oo. The magnitude of 
transverse magnetization \F+\ is proportional to 



tanh 
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Thus, the width of the domain wall is of the same order 
as the size of the spin vortex. 

The planar dark soliton is known to be dynamically 
unstable against "snake instability" [HI [Hi . This insta- 
bility arises for the wavelength longer than the critical 
wavelength, 
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After the distortion by the snake instability, the dark 
soliton transforms into vortex-antivortex pairs . This 
phenomenon has been observed in a nonlinear optical 
medium |2£} and in a two-component BEC [2]]. In the 
spinor BEC, an analogous instability causes formation of 
polar-core vortex-antivortex pairs. We will show in the 
next section that spin vortices are generated from domain 
walls by this mechanism. 



III. QUENCHING DYNAMICS OF THE 
TRAPPED SYSTEM 



A. The Berkeley experiment 

In this section, we will numerically simulate the mag- 
netization dynamics in a situation corresponding to the 
Berkeley experiment || and compare the results. We 
first briefly review the experiment. 

The trapping frequencies of the optical potential used 
in the experiment are given by (uo x1 uo yi uo z ) = 2tt x 
(56, 350, 4.3) Hz, and the system is effectively quasi-two 
dimensional (2D) in the x-z plane. The atoms prepared 
in the m = — 1 state are transferred to the m = state 
by rf field, and the magnetic field of B = 2 G is applied 
in the z direction. The number of spin-1 87 Rb atoms in 
the BEC is 2.1 x 10 6 with a peak density of n = 2.8 x 10 14 
cm -3 . These conditions give g/(|ci|n) ^ 28.4, and there- 
fore the prepared state is stable (see Fig. [I). The residual 
component in each of the m = ±1 states is less than 0.3%. 

The strength of the magnetic field is then suddenly 
decreased to 50 mG, which corresponds to q/(\ci\n) c± 
0.018. It follows from Fig. [I] that the m = state is 
no longer the ground state and spontaneous magnetiza- 
tion emerges. From t = 50 ms to 100 ms, the transverse 
(x-y direction) magnetization grows exponentially with 
the time constant of 15 ms. The spin vector varies in 
space and points in various random directions, and com- 
plicated spin textures can be observed. The polar-core 
vortices are identified in about one-third of the snapshots 
of the spin distribution, and sometimes several vortices 
coexist in a single sample. The correlation function for 
the transverse magnetization, 



G T (Sr) = Re 



/ drF_(r)F + (r + 5r) 



J drn(r)n(r + Sr) 



(24) 



oscillates in both x and z directions. The wavelength of 
the oscillation in the x direction is ~ 10 /am and that in 
the z direction is ~ 50 jam. The longitudinal magneti- 
zation Gl(0) shows no significant change within t < 300 
ms, where 



G L (Sr) 



J drF z (r)F z (r + 5r) 
J drn(r)n(r + Sr) 



(25) 



The quantities discussed in Sec.lTTIand observed in the 
experiment are summarized in Tab.[I] 



B. Numerical method 

We perform full 3D numerical calculations of the GP 
equation J7J) using the Crank-Nicolson scheme with a typ- 
ical grid size of c± 0.4 /am. In the numerical calculations, 
we ignore the linear Zeeman terms in the GP equation 
for the reason mentioned below Eq. J7J) . We first prepare 
the ground state in the m = state, i/jq = V^ni, by the 
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TABLE I: Length and time scales obtained by the mean- field 
theory and those observed in the Berkeley experiment [3. 
The last row of the right column is blank, since it has not yet 
been observed. 



mean-field theory 


experiment 


A mu in Eq. JTTJ 


wavelength of oscillation in Gr(Sr) 


~ 15.1 /mi 


~ 10 /im (in xj 




~ 50 /im (in z) 


T mu in Eq. (fT^|> 


time constant of Grr(0) 


~ 16 ms 


~ 15 ms 


£ s in Eq. flSJ 


size of vortex core 


~ 2.4 fim 


~ 3 /im 


Acr in Eq. 


length scale of vortex-antivortex pair 


~ 21 /xm 


? 



imaginary-time propagation of the GP equation (0 with 

= o. 

If the initial populations of the m = ±1 states are ex- 
actly zero, no spin-exchange dynamics follow within the 
GP equation. We therefore add small initial seeds to 
the m = ±1 states to trigger the spin-exchange dynam- 
ics. Possible physical origins of the initial seed include (i) 
residual atoms due to imperfections in the rf transfer, (ii) 
quantum fluctuations, and (iii) thermal components. In 
the experiment, although the residual fraction in each of 
the m = ±1 states is suppressed below 0.3% 8], this up- 
per limit corresponds to 6300 atoms, which is still large 
enough to significantly affect the subsequent dynamics. 
The spatial distribution of the residual atoms should re- 
flect that of the condensate = ^ini, while its phase and 
magnitude can fluctuate spatially due to experimental 
noise. Quantum fluctuations trigger the spin-exchange 
dynamics even if the initial populations in the m = ±1 
states are exactly zero, since the V^l^o term in the 
second quantized Hamiltonian transfers the m = pop- 
ulation to the m = ±1 states. The quantum fluctuations 
can be taken into account by random noises in the ini- 
tial state |22j |. The thermal component also triggers the 
growth of the m = ±1 states having phase fluctuations. 

In order to find an appropriate initial seed to repro- 
duce the experimental results and capture the essential 
mechanism that triggers the magnetization, we will ex- 
amine three kinds of initial seeds which reflect the shape 
of the condensate and various types of noise. We will see 
that the dynamics crucially depend on the nature of the 
initial seed. 



C. Initial seed without noise 

In the Berkeley experiment the atoms are first prepared 
in the m = — 1 state and then transferred to the m = 
state. We suppose that the transfer is imperfect and a 
small fraction is left in the m = — 1 state. Thus we 
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FIG. 2: Magnitude and direction of the spin at y = for the 
initial condition in Eq. (|^|> with e = 0.05. The upper panels 
show transverse spin density and the lower panels show 
its phase arg(i 7 +). The gray scale represents the density from 
to the peak value D p = n(r = 0) for the upper panels, and 
the phase from — tt to tt for the lower panels. The field of view 
of each panel is 400 x 25 /im. 



assume the initial state to be 

^i = 0, 

^ = V / T^^ini, (26) 
^-1 = ^ini, 

where e is a small constant. We take e = 0.05 and hence 
the initial population in the m = — 1 state is 0.25%, which 
is consistent with the experimental condition (< 0.3%). 

Figure El shows snapshots of the transverse spin density 
and its direction arg(F + ). The spin-density profile 
at t = 130 ms exhibits three magnetic domains, where 
the middle one magnetizes in the x direction and the two 
side ones in the — x direction. The width of each domain 
~ 8 /im roughly equals half of the most unstable wave- 
length pip. As time proceeds, new domains are formed 
at both ends of the cigar-shaped trap, whose alignment 
is perpendicular to the central domains (see the snapshot 
at t = 488 ms in Fig.EJ). Some of them become unstable 
and spin vortex-antivortex pairs are produced (enclosed 
by the square at t = 488 ms in Fig.EJ). Afterward, the two 
central domain walls become wavy (t = 590 ms), which is 
followed by the formation of spin vortex-antivortex pairs 
(enclosed by the square at t = 619 ms in Fig. EJ). Around 
the spin vortex the direction of the spin vector rotates 
by 2tt in the x-y plane, and the core is occupied by the 
m = component. This topological defect is therefore 
the polar-core vortex discussed in Sec. Ill C II Thus, the 
domain walls are first formed and the polar-core vortices 
then develop from the domain walls. 

We show that the dark-soliton picture discussed in 
Sec. Ill C 21 describes the dynamics shown in Fig. [3 very 
well. The domain walls generated in Fig. can be re- 
garded as dark solitons according to Eqs. ([T3|) and (|21|). 
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FIG. 4: Magnitude and direction of the spin for the initial 
condition given in Eq. (|7f)> . 



rection in Fig. 0(b) were not observed in the experiment, 
the initial seed in Eq. (|26|) does not correspond to that 
in the experiment. 



FIG. 3: (Color) (a) Time evolution of G T (0) (solid curve) 
and Gl(0) (dashed curve) and (b) Gt(x,z) for the initial 
condition given in Eq. with e 0.05. 

Fitting with Eq. (122)). we find the width of the do- 
main walls along the z axis at t = 130 ms is c± 2 /im, 
and ~ 3.3 /im for those along the x axis at £ > 420 ms. 
The latter is in good agreement with V2£, s — 3.4 /im and 
the former is somewhat smaller probably due to the in- 
fluence of the trapping potential. The snake instability 
manifests itself as the wavy domain walls at t = 590 ms 
with a wavelength ~ 19 /im, which is roughly equal to 
Eq. ([25|). 

Figure[3|(a) shows time evolution of Gt(0) and Gl(0), 
which indicate the degrees of magnetization in the x-y 
and z directions, respectively. The transverse magneti- 
zation Gt(0) exponentially grows in ~ 100 ms with a 
time constant ^16 ms, which is in good agreement with 
the experimental observation of 15 ms. The oscillation 
of Gt(0) at the frequency of ~ 8 Hz is seen from 100 ms 
to 600 ms, which was not observed in the experiment. 
The longitudinal magnetization Gl(0) remains small for 
t < 500 ms. The transverse correlation function Gr(Sr) 
is shown in Fig. [3| (b). The stripe pattern in the x di- 
rection emerges at t ~ 130 ms, which becomes unstable 
(^ 420 ms), transforming into the stripe in the z direc- 
tion. 

The initial growth of the magnetic domains staggered 
in the x direction is due to the anisotropy in the mo- 
mentum distribution in the initial seed (|26j) . The Fourier 
component of the most unstable wavelength (|11|) con- 
tained in the initial seed has the largest weight in the x 
direction. Since the domain walls staggered in the x di- 



D. Initial seed with white noise 

Next we examine the case of white noise as an initial 
seed in the m = ±1 states as 

^1 = M(h +^/4)/env, 

^o=A/> ini , (27) 

^-1 =A/*(/l +i/ 2 )/env, 

where AT is a normalization constant, f^s are 
random numbers obeying the normal distribution 
e -f /( 2cr ) J (y^^nr), and f env is an envelope function. The 
random number is chosen independently on each grid. 
We take a = 2 x 10 -4 , and the initial population in each 
of the m = ±1 states is ~ 0.23%. The envelope function 
J env is taken to be ^i n i. 

Figure 21 shows the distributions of the magnitude and 
direction of the spin for initial state (|27|). In contrast to 
Fig. the domain structure staggered in the z direction 
first emerges at t c± 170 ms. The size of the single domain 
is ^ 9 /im, which is roughly the same as the domain 
width in Fig. Oat t = 130 ms, reflecting the fact that the 
domain size is set by the most unstable length scale of the 
system. Some of the domain walls then develop into the 
polar-core vortices as shown in the squares in Fig. 21 and 
the system exhibits complicated spin dynamics similar to 
the experimental results described in Ref. p. 

Time evolution of the transverse Gt(0) and longitu- 
dinal Gl(0) squared magnetization is shown in Fig. 
(a). The time at which Gt(0) rises in Fig. [51(a) is later 
than that in Fig. [31 This is because the most unstable 
Fourier component in the initial seed is smaller in the 
present case due to the broad momentum distribution 
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FIG. 5: (Color) (a) Time evolution of G T (0) (solid curve) 
and Gl(0) (dashed curve) and (b) Gt(x,z) for the initial 
condition given in Eq. ((27)) . 



of the white noise. Figure [51 (b) shows the transverse 
correlation function Gt(Sv). The clear stripe pattern at 
t = 204 ms indicates the establishment of the long-range 
correlation over > 200 /mi. The stripe pattern is then 
distorted and its width becomes broader for t > 600 ms. 

We note that the qualitative behaviors in Figs. ^ and 
are relatively insensitive to the envelope function / env 
of the initial seed in Eq. ((27)). We have used / env = ^i n i 
based on the assumption that the noise reflects the shape 
of the condensate. However, we find that the results re- 
main qualitatively the same even when the envelope func- 
tion is not multiplied. We have also confirmed that the 
results shown above are insensitive to grid size, despite 
the fact that a random number is initially assigned to 
each grid and the momentum distribution of the initial 
noise depends on the grid size. 



E. Initial seed with colored noise 

The initial states (|2(j[) and ((27)1 examined in Sees. 1111 
and IIII Dl have been inadequate to reproduce the exper- 
imental results. In the experiments, the transverse cor- 
relation function Gr(Sr) oscillates both in the x and z 
directions and the wavelength of the oscillation is larger 
in the z direction than that in the x direction. For the 
initial seed proportional to ^i n i given in Eq. ((25)). the do- 
mains staggered in the short axis first grow, and for the 
white noise in Eq. ((27)) . the domains are staggered in the 
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FIG. 6: Magnitude and direction of the spin for the initial 
condition given in Eq. (|28f) with A cu toff = 60 /im. 



long axis. Therefore, we expect that the experimental 
result can be reproduced using the initial seed that com- 
bines Eqs. ((26)) and ((27)) . i.e., anisotropy and randomness. 
Since the white noise in Eq. ((27)) induces the growth of 
the domains along the x direction (Fig. ^ t = 170 ms) 
and the wavelength of the oscillation of Gr(Sr)^ 18 
/im (Fig. [SJ t = 204 ms), is shorter than that in the 
experiment, ~ 50 /im, we cut off the short- wavelength 
components from the white noise. The initial state is 
thus given by 



4>- 



: AA/'ini, 

= -A/" /cutoff /env i 



(28) 



where we produce the noise function /cutoff from the 
white noise by eliminating the Fourier components whose 
wavelengths are shorter than A cuto ff. The envelope func- 
tion / env is taken to be ^hii- 

The results for A cuto ff = 60 jam are shown in Figs. Eland 
[3 The behavior of Gt(Sv) shown in Fig. 0(b) is similar 
to that observed in the experiment 8] in that the wave- 
length of the oscillation in the z direction is larger than 
that in the x direction. This anisotropy originates from 
the envelope function / env = of the initial seed in 
Eq. ((28)). The momentum distribution in the x direction 
of /env is broad and the most unstable wavelength (|llj) 
grows in this direction, while in the z direction the typ- 
ical wavelength in Gr(Sr) is determined by the cutoff 
wavelength A cuto ff of the noise. 

We note that the qualitative behaviors shown in Fig.0 
(b) are independent of the details of the envelope function 
/ env in the initial seed. The similar behaviors can be 
obtained as long as the size of the envelope function is 
much larger than the unstable wavelength A mu in the z 
direction and is comparable to A mu in the x direction. 
In fact, when we use / env = Vfni? qualitatively similar 
results are obtained. 
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FIG. 7: (Color) (a) Time evolution of G T (0) (solid curve) 
and Gl(0) (dashed curve) and (b) Gt(x,z) for the initial 
condition given in Eq. (|28f> with A cuto ff = 60 /im. 
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FIG. 8: Time evolution of the system, in which the magnetic 
field is decreased from 530 mG to 50 mG during the first 300 
ms so that q linearly decreases. The initial condition is the 
same as m Figs. El and 
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The behavior of the local transverse magnetization 
Gt(0) in Fig. 0(a) is also in good agreement with the ex- 
perimental result. It exponentially increases from t = 50 
ms to 100 ms with a time constant ~ 15 ms. After 
t = 100 ms, Gt(0) gradually increases from 0.4 to 0.6 
until t = 300 ms, which also captures the basic charac- 
teristics of the experimental result. 

As in Figs. and 21 Fig. El shows that the magnetic 
domains are first formed (t = 102 ms) followed by the 
development of some of the domain walls into the polar- 
core vortices. Thus, this process of spin-vortex formation 
appears rather universal. We also found that the polar- 
core vortices drift in and out of the condensate in the 
dynamics. 



F. Dependence of the dynamics on quench time 
and magnetic field 

We have so far considered the case of sudden quench, 
i.e., the magnetic field being suddenly reduced to 50 mG 
at t = 0. If the time scale of the quench is longer than 
^ 100 ms (time scale in which Gt(0) rises in Fig.fTJ), the 
excitations are expected to be suppressed because of the 
adiabatic theorem. 

Figure [HI shows the time evolution for the slow quench, 
where the initial magnetic field of 530 mG is reduced to 
50 mG during 300 ms so that q linearly decreases. We 
find from Fig. 0(a) that the spin state has nearly a single- 
domain structure and no spin vortices are created. The 
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FIG. 9: Ramp time dependence of the number of spin vortices 
at t = 200 ms. The initial state is given in Eq. (|28f) with 
Acutoff = 60 /im and the magnetic field is decreased from 530 
mG to 50 mG so that q is linearly ramped down during tramp • 
The plots and error bars represent the average and standard 
deviation with respect to ten runs of simulations for different 
random numbers to generate the noise. 



transverse Gt(0) and longitudinal Gl(0) components of 
the squared magnetization are shown in Fig. EI (t>). A 
large value of Gt(0) ~ 0.9 is due to the absence of the 
spatial spin structure. 

Figure El shows the dependence of the number of spin 
vortices at t = 200 ms on ramp time. The number of spin 
vortices decreases with an increase in the ramp time, and 
the time scale beyond which the spin vortices no longer 
emerge is given by ~ 100 ms. We note that this time 
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FIG. 10: Time evolution of the system, in which the magnetic 
field is suddenly decreased to 400 mG at t = 0. The initial 
condition is the same as in Figs. El and □ 



scale coincides with /i/(|ci|n), indicating that the energy 
scale for the creation of a spin vortex is given by ~ \c\\n. 
As in the experiment, the number of vortices fluctuates 
from run to run. 

The result seen in Fig. El reminds us of the Kibble- 
Zurek mechanism 0, 0] of the vortex creation in the 
quenched system, which also depends on the quench 
time. However, the relationship between the present 
phenomenon and the Kibble- Zurek mechanism is not 
straightforward. In the Kibble-Zurek mechanism, each 
domain is assumed to be created with an independent 
phase. In the present system, however, the magnetization 
at each position cannot be independent of that of other 
positions, since there is the restriction on the change in 
the total spin, which is clearly seen in the long-range cor- 
relation in Gr(Sr). Gaining a deeper understanding of 
the relationship between the spin-vortex formation and 
the Kibble-Zurek mechanism constitutes an interesting 
and challenging problem which merits further investiga- 
tion. 

The spin-vortex formation is also suppressed if the 
quench of the magnetic field is made just below the crit- 
ical strength of the magnetic field, which is given by 
q = 2|ci|n (see Fig. and corresponds to B = 530 
mG for the peak density. Figure E3 shows the spin 
dynamics for B = 400 mG. From Eq. (|T0j) . the most 
unstable wavelength for this magnetic field is A = oo, 
and the system is more unstable against excitations with 
larger wavelengths. In Fig. [HJ (a), therefore, only the 
long-wavelength modes are excited and no spin vortex 
is formed. In Fig. (b), the transverse squared mag- 
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FIG. 11: Time evolution of the transverse magnetization 
= | J drF + \ for the dynamics in Fig. |3 (solid curve), 
Fig. 2] (dashed curve), and Fig. |S] (dotted curve). 



netization Gt(0) saturates around 0.2, since \F+\ of the 
ground state is small near the phase boundary between 
the broken- axisymmetry and polar phases. 



G. Total transverse magnetization 



We define the total magnetization as 



(F) 



drF, 



(29) 



where F is the spin density defined in Eq. ©. If the 
quadratic Zeeman effect is absent, the z component (F z ) 
and the transverse component are conserved, and 

the vector ((F x ),(F y )) rotates in the x-y plane at the 
Larmor frequency determined by the linear Zeeman en- 
ergy. In the presence of the quadratic Zeeman effect, not 
only (F) rotates in the x-y plane but also \(F+}\ changes 
with time. 

Figure ITT1 shows time evolution of \(F+}\ for the initial 
conditions given by Eqs. (|2T)h and J2EJ), where the 

system is quenched by a decrease in the magnetic field 
to 50 mG at t = as in Figs. 1-171 For the initial con- 
dition (J2HJ), \{F+)\ has a large initial value ^ 0.07, and 
reaches a maximum value ~ 0.19. For the initial con- 
dition (|28j). | monotonically increases to twice the 
initial value during 600 ms. Thus the quadratic Zeeman 
effect generates the transverse component in the total 
magnetization. 

The transverse magnetization may also be 

changed by the dipole-dipole interaction, since it cou- 
ples the spin angular momentum with the orbital angular 
momentum |23|, |24| . The effect of the dipole-dipole in- 
teraction on the transverse magnetization merits further 
study. 
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IV. CONCLUSIONS 

We have studied the spontaneous magnetization and 
spin-texture formation of a spin-1 87 Rb BEC, where the 
initial state is the m = stationary state ^i n i plus a small 
seed in the m = ±1 states. 

We have reproduced the polar-core spin- vortex forma- 
tion observed in the experiment, as shown in Figs. [2j 
and El Typically, the spin vortex is formed in two steps. 
The magnetic domains are first formed, and then the 
domain walls transform into the spin vortex- ant ivortex 
pairs. This process of vortex formation appears to be 
universal regardless of various initial conditions. 

We have examined three kinds of initial seeds: the one 
proportional to ipi n {, white noise, and colored noise. The 
first two seeds produce the domain structures staggered 
in the short and long axes, respectively (Figs. The 
magnetization developed from the third seed has both 
characteristics of the first two seeds, and the correlation 
function oscillates in both the long and the short axes 
(Fig. 13(a)), in qualitative agreement with the Berkeley 
experiment Q . This is due to the fact that the third seed 
has a broad momentum distribution in the short axis, 
originating from the shape of ip- m i, and long- wavelength 
fluctuations in the long axis. From these results, we 
can conclude that the anisotropy and colored noise in 
the initial seed are important to account for the exper- 
iment. The time evolution of the transverse magnetiza- 
tion Gt(0) is also in close agreement with the experimen- 
tal result (Fig. 13(b)). 

The number of spin vortices created in the magnetiza- 



tion depends on how fast the magnetic field is quenched. 
When the magnetic field is decreased slowly, the num- 
ber of nucleated spin vortices decreases (Figs. |EI and EJ). 
We have also shown that the number of spin vortices de- 
creases for the magnetic field close to the critical value 
(Fig. EH . 

The transverse component of the total magnetization 
\(F+)\ is changed by the quadratic Zeeman effect, and 
can exceed twice the initial value (Fig. ITT]). 

We have pointed out a close analogy between the topo- 
logical defects in the present system and those in a scalar 
BEC, and that the creation of spin vortex- ant ivortex 
pairs from the domain walls is related to the instability 
in the planar dark solitons. It is of interest to investigate 
if the counterparts of vortex lattices, multiply-quantized 
vortices, and gray solitons are generated in a ferromag- 
netic spinor BEC. 
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